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We demonstrate which way slowing down ultrarelativistic shocks create GRB lags. Reflection process produces positive 
lags and Cracow acceleration produces negative lags. We describe the process of Cracow acceleration and present two 
ways in which the seed particles are injected into the upstream plasma. Strong decelerating shocks in the presence of 
the acceleration processes account for the observed hard energy spectra of accelerated electrons with spectral indices 
smaller than the value 2.2. We claim that during the strong deceleration stage the rise-fraction of seed particles is formed 
upstream of the shock. The rise-fraction feeds the Cracow surge and normal seed particles feed the reflection surge. We 
present the model of the microphysics of relativistic plasma that aims at explaining the required disturbances of the 
movement of particles upstream of the shock that allow for Cracow acceleration. We show that Cracow acceleration 
can produce UHECRs. 

Key words, acceleration of particles — shock waves — cosmic rays — radiation mechanisms: non-thermal 



1. Introduction 



iFermil (|1949T ) proposed a mechanism being supposed to ex- 
plain cosmic rays acceleration. The idea is analogous to the 
acceleration of a tennis ball. Interstellar clouds are tennis 
rackets and a particle is the tennis ball. There is a small 
difference here, tennis rackets rarely decrease the energy 
of the ball. Cosmic rays are accelerated in this mechanism 
according to second order Fermi acceleration method. It 
should be called Fermi acceleration. 

The concept that shock waves could acc elerate parti- 
cles appeared slowly. It was foreshadowed bv iHovlel (i960) 
who postulated that shocks could efficiently accelera t e par - 
ticl es but w i thout specif ying a mechanism. iParkerl ()1958[ ) 
and iHudsonl (|1965L Il967li attempted to obtain such mecha- 
nism based on pai rs of converging shocks and, most notably, 
iSchatzmanl (|1963f ) constructed a theory based on perpendic- 
ular shocks where particles keep bumping into a hydromag- 
netic shock front and increase their energy in the Fermi-like 
way. 

The real ac celeration mecha n ism was descr i bed in four 
semina l papers, iKrvmskvl (ll977l).lAxford et al.l ljl977f ).lBel3 
(jl978allbf ) and iBlandford fc Qstrikerl (| 19781 ). Since many 
persons aspired to the name of the discoverer, therefore 
this mechanism was called diffusive shock acceleration. It 
is often called Fermi acceleration, but it is incorrect, since 
Fermi does not have anything to do with particle acceler- 
ation at shoc ks. The name stuck to this mechanism (see 
lHoshinoll2008l page 940) and we use it in this paper. 

One should, by the way, recall nomenclatures that mis- 
lead uninformed readers about a misunderstanding which 
persists in the literature. First order Fermi acceleration and 
second order Fermi acceleration are not a real mechanism. 
It is the method which produces the Fermi statistical de- 
pendence in which the proportional gain of energy is the 
same for all particles. There are two or more media in the 
method that are moving with the average speed equal to 



V. The speed of a particle is much larger that V and is 
equal to c. The particle bounces between the media and 
gains energy. If the direction of V is always opposite to the 
direction of c in each collision then the average energy gain 
of the particle is proportional to V/ c and that is first order 
Fermi acceleration. If the directions are random then the 
average energy gain is proportional to {V/c) 2 and that is 
second order Fermi acceleration. The both methods appear 
in separate phenomena and there is no need of applying a 
shared name for them. It could sometimes lead to ambi- 
guity. In rare cases when it is necessar y one can us e the 
phrase 'all Fermi acceleration processes' (|Drurv| |1983. page 
987) . First order Fermi acceleration assumes that V <C c. It 
implies that the distribution of accelerating particles in the 
medium becomes nearly isotropic so that the particles could 
reach the opposite medium. First order Fermi acceleration 
predicts that the energy of accelerating particles increases 
~ 7 2 times in each cycle upstream-downstream-upstream 
if the mechanism is applied to relativistic shocks. 

Some authors are trying to name the idea of the particle 
energy increase as a result of particles bouncing back and 
forth across the shock wave as Fermi acceleration, but it is 
incorrect. It is Hoyle-Schatzman idea. 

In Fermi acceleration particles can cross the shock re- 
peatedly thanks to magnetic field fluctuations downstream 
of the shock and subluminal shock geometry. Magnetic field 
fluctuations upstream of the shock do not play significant 
role here, since the shock will always catch up with the 
particle that is wandering there. This mechan ism does not 
accelerate pa rticles in superluminal shocks (|Belll . Il978at 
iDrurvl Il983|) and therefore is limited to non-relativistic 
and poorly relativistic shocks, since relativistic shocks are 
generically superluminal. 

For a lot of years an unresolved riddle existed in astron- 
omy. Observations have shown that in areas where relativis- 
tic shocks are present particles are accelerating effectively, 
but nobody was able to find the mechanism. 
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Fig. 1. Fermi acceleration arises due to inclination and fluc- 
tuations of the magnetic field downstream of the shock. 
Cracow acceleration arises due to fluctuations of the field 
upstream of the shock. 



The acceleration mechanism that operates at 
ultrarelativistic sho ck fro nts was discovered by 
iBednarz fc Ostrowskil ([19981 ). From now on we will 
call it Cracow acceleration since it was discovered in 
Cracow. This imp ortant discovery was p r esente d earlier at 
three conferences, Bednarz & Ostrowski (1997a. b.c). 

IBednarz fc Ostrowskil (|1998|) showed that ultrarelativis- 
tic shocks accelerate particles and accelerate them more ef- 
fectively if their Lorentz factor grows. It results from the 
place and the way the particles are getting the chance of 
the return to the shock. In Fermi acceleration they return 
to the shock because of conditions downstream of the shock 
and in Cracow acceleration because of magnetic field fluc- 
tuations upstream of the shock (Figure [T|) . 

There were a few proposals before, but th eir mech- 
anisms are ineffective. iBeeelman fc Kirkl (|l990t l proposed 
shock-drift acceleration at relativistic shocks to operate 
at hot spots o f extr agalactic radio sources. Afterwards, 
iHoshino et al.l (|1992T ) described a process of shock ac- 
celeration of positrons to non-thermal distributions de- 
voted to account for the synchrotron radiation observed 
in the Crab Nebula and h ot spots. Other proposals ar e 
magnetic reconnection (e.g. iRomanova fc Lovelace! [T 992') . 
shock surfing acceleration ( e.g. ISagdeev fc Shapiro! 119731 
Katsouleas fc Dawson! Il983l ) and MHD turbulence (e.g. 



Bieber et al!I Tl994). 



Reflection process is connected with relativistic shocks. 
It was known that a particle increases its ener gy 7 2 
times in the cycle upstream-downstream-upstream. Vietril 
(|l995l) proposed repeated crossings at an ultrarelativistic 
shock to occur and argued that it would lead to ^ en 



ergy gain per crossing cycle. IBednarz fc Ostrowsk i (1998) 
found that because of strong anisotropy of cosmic rays 
upstream of the shock the gain equals to ~ 2, what was 
confirmed bv lGallant fc Achterberd (|1999D . However, flow- 
ing in particles have isotropic distribution and gain j 2 
times energy in the fir s t cycl e and it was called reflection 
(jBednarz fc Ostrowskil Il999h . In this paper we will show 
that reflection process is an important element of the cos- 
mic ray acceleration. 

After IBednarz fc Ostrowskil (|l998l ) they have been try- 
ing to include the idea of cosmic ray anisotropy up- 
strea m of the shock, b u t in the Fermi ac c elerat io n con- 
text, iKirk et al.l ([20001 ): lAchterberg et al.l (|200ll) : IVietril 



( 200 3h ; iLemoine fc Pelletien ([20031 ) - disordered shocks and 
Ell ison fc Double! (120021 ) - parallel shocks. Fermi acceler- 
ation takes place at relativistic shocks in two cases not- 
appearing in the outer space, at practically parallel shocks 
(subluminal geometry) and at shocks with completely dis- 
ordered magnetic field downstream of the shock (we named 
them disordered shocks) . One would call these two theoreti- 
cal cases Peacock acceleration ([PeacockM 19811 ). The Weibel 
instability is preventing disordered and parallel shocks from 
being form ed what observations of polarised em ission are 
co nfirming (ISteele et all l2009t lG5tz et alj, |2009T>. 

Niemiec fc Ostrowskil (|2006l ). iNiemiec et all ([2006D . 



ILemoine et al.l (|2006l) were trying to examine the accelera- 
tion process and used 'realistic' magnetic field fluctuations 
upstream of the shock. They have failed to get e f fective 
acceleration at superluminal shocks. iPelletier et al.1 ([2009) 
have been discussed it analytically and have been found 
to agree with the numerical results. They have concluded 
that the acceleration takes place if the fluctuations ampli- 
tude SB is much larger than the homogeneous magnetic 
field Bq upstream of the shock what produces disordered 
magnetic field downstream of the shock, but this is the case 
of Peacock acceleration. Their lack of effective acceleration 
results from this that they have applied normal conditions 
(types of fluctuation spectra) to the phenomenon which is 
completely unknown. 

Cracow acceleration needs a type of the fluctuation 
spectrum where almost entire fluctuation energy is gath- 
ered in very intense waves about a small length. The length 
can be so small as the size of the atom. We do not belive 
in such fluctuations and propose an other solution to the 
problem in Section [5] 

We expect that astronomical sites where particles are 
accelerated in Cracow acceleration are relativistic shocks 
with the Lorentz factor much larger than 1 and the helical 
(perpendicular to the shock normal) mean magnetic field 
prod uced by the Weibel instability (jMedvedev fc LoeR 
1999). In GRBs the process does occur not only in prompt 
gamm a-ray pulses, but als o in X-ray flares dMargutti et all 
2010) and m precursors ([Burlon et all I2009D where" very 
similar spectral properties are observed. Long and short 
GRBs are also sho wing the same acceleration process 
dGuiriec et all 120101) . 

In pulsar wind nebulae (PWNe) where (contrary to 
GRBs) the plasma downstream of the shock is approxi- 
mately at rest with respect to ISM, not only lept o ns bu t 
probably protons are being accelerated (|Li et all [2010). 
There not only the same mechanism of the acceleration 
is present, but probably similar physical conditions, since 
the assumed particle injection spectrum in the form of a 
broken-power law turns out t o have spectral break a t an 
similar energy for all sources ([Bucciantini et ail . 1201 If) . 

We think that particles are accelerated in the mech- 
anism of Cracow acceleration in gamma-ray binaries 
(jCerutti et all 120091 ). in X-ray binaries where relativistic 
jets are formed and in jets thrown away by active galactic 
nuclei. An example of X-ray binaries is Cygnus X-3 where 
electrons gain energy at the place where the jet is recol- 
limated by the ste llar wind pressure and forms a shock 
(|Dubus et all 120101 ) . 

In the test particle limit the energy spectral index of 
accelerated particles is only dependent on the compression 
ratio of the p l asma through the shock by /3 = (R+2) / (R— 1) 
(|Bell Il978at iDrurvl Il983h . The formula is not valid for 
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relativistic shocks but it is good enough for the estimation 
how j3 increases when the compression falls down. 

The magnetization parameter, a, is equal to the up- 
stre am Poynting flux relative to the total mass-energ y 



stre am .royntmg f lux rel ative to trie total mas s-energy 
flux (lAppl fc Camenzindl 119881: Kennel fc Coron"itlll984a|K 



iKennel fe CoronitiM 1984a|) presented the derivation of the 
Rankine-Hugoniot relations for perpendicular shocks. They 
found that R = 3(1 - 4a) for a < 0.01 and R = 1 + l/(2cr) 
for a > 10. We estimate that Cracow acceleration needs 
a < 0.1 — 0.01 to accelerate efficiently if one applies the 
usual plasma conditions, but we devise the model of loop 
plasma which could require large value of a and standard 
co mpression R = 3. 

iKirk fc Ski i] (|2003h have found that the Poynting 

flux can be dissipated before the pulsar wind reaches the 
inner edge of the Crab Nebula. This is in accordance with 
the val ue of a = 0.003 near the t ermination shock of the 
nebula (jKennel fc Coronltll . 1984a. b) which is far more suf- 
ficient than Cracow acceleration requires. 

The context of the paper is organized as follows. In 
Section [2] we present the model of the microphysics of rela- 
tivistic plasma as well as the way the simulations of the rela- 
tivistic plasma should be conducted. In Section[3]the mech- 
anism of Cracow acceleration and in Section0]thc numerical 
simulations are described. In Section [5] we give thought to 
the influence of the deceleration of the shock on the energy 
spectra of accelerated particles. Results of the simulations 
are presented in Section [6l there are many tables and fig- 
ures. Section [7] is devoted to the problem in what condi- 
tions Fermi acceleration is able to produce negative lags. 
In Section [5] we postulate the formation of the rise- fraction 
of particles upstream of the shock at the strong deceler- 
ation stage. The rise-fraction feeds the Cracow surge and 
normal seed particles feed the reflection surge. In Section [3] 
we show that Cracow acceleration can produce UHECRs. 
The seed particles problem is discussed in Section [TUJ We 
summarize and discuss our paper in Section QTJ 



2. The microphysics of relativistic flows 

Relativistic plasma is an obscure phenomenon. Magnetic 
waves are too weak to disturb the particle movement so 
strong as Cracow acceleration requires. There must be some 
collective processes that build magnetic and electric fields 
on a micro-scale. Before we introduce the model of loop 
plasma let us present how small magnetic helical tubes 
could build the homogeneous magnetic field and the re- 
quired disturbances. Let us consider two idealised cases. 

Let us assume that the Weibel instability generates on a 
micro-scale small magnetic helical tubes. The axes of sym- 
metry of these tubes are parallel to the global speed of the 
plasma. The tubes have different diameters and are scat- 
tered randomly across the space. Let us assume that the 
direction of the velocity of accelerated particles can diverge 
only a bit from the direction of axes of the tubes. This is 
the case of the accelerated particle wandering upstream of 
the relativistic shock. When the particle crosses the trans- 
verse section of a tube then the direction of its movement 
is changing slightly. The change is random and its value 
and direction depend on the location of the intersection 
point where the particle crosses the transverse section of 
the tube. These crossings correspond to scatterings in our 
simulations. 



In addition, we assume that vectors of the magnetic 
field of small tubes are adding up and form one magnetic 
helical tube in the macro scale (the big tube). The particle 
wandering upstream of the shock can see the big tube as 
the homogeneous magnetic field perpendicular to the shock 
normal. Now, we have two idealised cases, the random field 
and the homogeneous field, that are one physical reality. 
The value of the magnetic field of a small tube must be on 
average much larger than the value of the magnetic field of 
the big tube so that the addition of vectors is correct. 

Now, our model of simulations fits into the theoretical 
model. In the theoretical model particles interact with the 
strong random magnetic field only and the homogeneous 
magnetic field is formed by the composition and dynamics 
of small tubes. Each scattering takes up a very short time 
span, equal to zero in the simulations, and corresponds to 
the change of the direction of the particle velocity caused 
by some group of small tubes. The wandering in the ho- 
mogeneous magnetic field between scatterings takes up the 
entire amount of the time (the big tube) . Our model of sim- 
ulations enables any shortening of the time between scat- 
terings and keeping the same pattern of the magnetic field 
fluctuations simultaneously. 

One should remember that current numerical simula- 
tions and analytical models are useless in the topic of rel- 
ativistic plasma since they do not reach the required ac- 
curacy. We will explain this opinion with the use of two 
ex amples of PIC sim ulations. 

Spitkovskv (2008) simulated relativistic shocks in 2D 
and used an artificial wall that had produced an artificial 
returning stream. Collective processes in plasma are sensi- 
tive to a flap of a bird's wing and the stream is a hurricane 
which can do everything. The author has written that the 
shock moves to the right at 0.47c. We do not know how to 
transforme the velocity to 3D, but if it has the same value 
in 3D then the compression is equal to 2. Such decline in 
the value of the compression from 3 to 2 means that the 
asymptotic value of the energy spectral index /3q grows from 
2.2 to ~ 4. The author has detected particle acceleration 
and got the index equal to 2.4. In our opinion the artifi- 
cial returning stream accelerates particles and reduces the 
co mpression. 

Nishikaw a et alj (|2006f) performed simulations with ~ 
3.8 • 10 s particles, but we have counted ~ 2.5 • 10 s what is 
unimportant here. The important thing is that 1.25 • 10 8 
hydrogen atoms have mass 2.1 ■ 10~ 16 g, 1.25 • 10 s electron- 
positron pairs have mass 1.1 • 10 _19 <7, a flu virus has mass 
7 • 10 _13 5 and typical plasma flows of GRBs have rest mass 
~ 10 27 <7 if we apply canonical energy release of ~ 10 51 erg, 
the Lorentz factor of the plasma equal to 1000 and neglect 
the Poynting flux. 

The authors have scaled their simulations to density of 
the interstellar medium in order to get the length of the 
simulation box equal to 4 • 10 cm and have not scaled to 
density of the tail of a comet, for example, which is 10 5 
times higher. The authors have compared the length of 
the simulation cell to the electron skin depth in order to 
show that the shock fits in the simulation box, but it is 
untrue. The thickness of th e shock is co mparable to the 
gyroradius of a thermal ion (|Drurvl . 119831 ) and it compares 
unfavourably with the size of their simulation box. The au- 
thors have not reached dynamic equilibrium so that they 
do not know if their results are stable. The authors have 
written, '.. the 'flat' (thick) jet fills the computational do- 
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Fig. 2. The model of loop plasma. The six solid line ellipses 
represent the upper layer and the six dashed line ellipses the 
bottom layer. The detailed description is in text. 



main in the transverse directions (infinite width). Thus, we 
are simulating a small section of a relativistic shock infi- 
nite in the transverse direction.' and that is untrue again. 
They have applied some boundary conditions and named 
them 'thick jet' or 'infinite width'. All boundary conditions 
are artificial and can produce anything. Further below, we 
propose how to manage with the problem. 

We have shown that the current PIC simulations of 
relativistic plasma are pointless. They have provided two 
intuitive hints only. First, relativist i c pla s ma could pro- 
duce seed particles (jNishikawa et all 120031: iHededal et all 
12004) and second, the pla sma could form current filaments 
(jMedvedev fc LoebL fl999l) . We believe that the task of cre- 
ation of relativistic plasma in PIC simulations from scratch 
is outside present capabilities of mankind. We come up with 
an idea and show the way the PIC simulations should ap- 
proach the problem of relativistic plasma. 

At first, we advance a hypothesis on the microphysics of 
relativistic plasma. The hypothesis aims at explaining the 
required disturbances of the movement of particles UP that 
allow for Cracow acceleration. Now, we restrict our con- 
siderations to the electron-positron plasma. Our model is 
presented in Figure [2j The Weibel instability forms electric 
current filaments in the shape of an oval which we substi- 
tute for an ellipse. Sometimes, we are calling this oval the 
loop. The sense of the electric current density j is the same 
for all loops. Each loop contains identical electric charges, 
electrons or positrons, which are indicated in Figure [2] by 
e~ , e + , respectively. The particles are arranged evenly along 
the loop in this way that the centrifugal force and the elec- 
tric force are balanced by the Lorentz force. The isolated 
loop in the shape of a circle is in ustable equilibrium and 
the environment flattens it into the oval. We have no idea 
which oval is more stable, the one elongated parallel or per- 
pendicular to the shock normal, but PIC simulations sug- 
gest that the ovals elongated parallel to the shock normal 
make relativistic plasma because current filaments tend to 
be parallel to the direction of the flow. 

We think that the system presented in Figure [2] is stable 
and the Weibel instability rebuilds continuously the ovals, 
but the PIC simulations are required to prove it. The sys- 
tem is certainly not so perfect as presented in Figure[2] The 
ovals have different diameters, they oscillate, have global 
movements (the instantaneous velocity of a loop is v e ). The 



ovals are being distorted in space. There is a fraction of par- 
ticles that do not belong to any oval (free particles). 

The ovals form the homogeneous magnetic field - B. 
They are divided into layers which are parallel to the shock 
normal. Each layer is divided into rows which are perpen- 
dicular to the shock normal. Ovals along a row are arranged 
alternately, the electron oval follows the positron oval and 
vice versa. We adopt the rule that the ovals approach the 
geometry that minimize the magnetic field opposite to B. 
We have not carried out any calculations, but we think that 
two adjacent rows adhere to the shared adhesive line. We 
indicate lengths of semi-axes of the ellipse by a and b, the 
distance between two consecutive ellipses along the row by 
d. Each layer is translated in regard to the symetrical po- 
sition of the adjacent layer for the length of a parallel to 
the shock normal and d/2 perpendicular to the shock nor- 
mal. The distance between two adjacent layers depends on 
values of a, b and d. 

In our model not only the magnetic helical tubes (Figure 
[5]) but also separation of electric charges disturb cosmic 
rays. It is easy to prove whether such disturbances en- 
able cosmic rays to accelerate in Cracow acceleration. One 
should track the trajectory of high-energy particles up- 
stream of the shock (Figure [5]) with their initial momenta 
nearly parallel to the shock normal. The change of the di- 
rection of the momenta must be small while the tracking. 
The momenta must be all the closer to the shock normal 
the Lorentz factor of the shock is larger. Now, one derives 
the value of Qq from the tracking and puts it into the for- 
mula [TJ If the derived energy spectral index is close to the 
asymptotic energy spectral index then the system acceler- 
ates particles in Cracow acceleration. 

The proton-electron plasma is more difficult to analyse. 
We think that it consists of proton loops and electron loops. 
Diameters of the electron ovals are much smaller then pro- 
ton ovals. The ovals also approach the geometry that mini- 
mize the magnetic field opposite to the homogeneous mag- 
netic field. Electron ovals tend to be close to the electric 
current of proton loops. 

We name the described plasma loop plasma. The dis- 
tribution of particle velocities in loop plasma is not the 
Maxwellian one so that the temperature is not defined in 
the usual meaning. It is hard to say what the downstream 
plasma will look like if loop plasma comes through the 
shock. However, it is not important for Cracow accelera- 
tion which will work anyway. 

The PIC simulations are useless unless one is able to 
detect dynamic equilibrium of plasma. Creating relativis- 
tic plasma from scratch is outside the available computing 
power at the moment. The PIC simulations must hunt for 
dynamic equilibrium of plasma rather than try and create 
it from scratch. Below, we describe how the simulations 
should be conducted. We will demonstrate how to catch 
dynamic equilibrium of electron-positron loop plasma. At 
first, we take into account the upstream loop plasma only. 
The shock is far enough and does not affect the system. 

One must use many processors at the same time (par- 
allel programming). The plasma is divided into boxes. The 
boxes are within the space of a cuboid (the central cuboid). 
Each electric current loop is within one loop box. Inside the 
loop box there are free particles also. Loop boxes must not 
overlap so that their linear dimensions are changing. If nec- 
essary, a few loops can be inside one box. The loop boxes 
can merge, split, disappear and be formed. Additionally, 
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the plasma is divided into free boxes that contain the re- 
maining free particles inside empty spaces. All free boxes 
are the same and fill up the cuboid entirely, they do not 
change. 

To save the computing power, one processor handles 
many loop boxes or free boxes during one computational 
step, but not at the same moment. The number of all boxes 
should exceed the number of processors ten times for ex- 
ample. The electromagnetic force inside a box is calculated 
exactly, the force from nearby boxes is also calculated ex- 
actly, but the force from any other box is taken from a 
pattern database. 

We replace the boundary conditions with a pattern 
database of the electromagnetic field of pattern cuboids. 
The pattern cuboids surround the central cuboid up to 
a distance. There are no empty spaces inside the system. 
The distance must be large enough to not change results 
of simulations if it increases. If necessary, global electric 
and magnetic fields that represent the electromagnetic field 
of plasma outside the distance is added. Sizes of pattern 
cuboids can be different and the results can not depend on 
the sizes. 

The biggest challenge is creating the pattern database 
and the efficient use of it. Each pattern in the database must 
have many items of accuracy. Let us say, we have 15 items 
of accuracy and number represents the lowest accuracy. 
Results of simulations must be the same if we apply items 
numbered from a certain number to number 14. 

The pattern database changes as simulations bring 
closer to the dynamic equilibrium target. The size of the 
central cuboid must be large enough to detect statistically 
independent pieces of it. Pattern cuboids will be formed 
by splitting and merging the central cuboid in order to get 
different sizes of them. 

Plasma inside the central cuboid moves with a velocity 
in the observer's frame. The plasma is supplied with boxes 
at one end of the central cuboid. The boxes disappear at 
the opposite end. The entering boxes are changing. They 
are replaced with boxes similar to those ones nearby the 
opposite end, but not too close to the boundary of the cen- 
tral cuboid. The simulations must run until all boxes inside 
the central cuboid are statistically homogeneous, pattern 
cuboids are similar to the central cuboid and the conditions 
do not change. It is the dynamic equilibrium of plasma. 

The next stage is to catch equilibrium of the plasma 
around the shock. It is easy because we know what the en- 
tering boxes look like. The central cuboid must contain the 
upstream plasma, the shock and the downstream plasma. 
Initial boxes of the shock and downstream plasma can be 
any. Simulations are conducted in the shock rest frame. 

The equation of motion of the shock must be known 
from astronomical observations in order to catch equilib- 
rium of the slowing down shock. Simulations should be con- 
ducted in the shock rest frame. The entering boxes are sup- 
plied according to the equation of motion. The simulations 
must start with different starting points that are between 
the initial and the final Lorentz factor of the shock. At the 
end, the plasma must change in the same way indepen- 
dently of the choice of the starting point. 

We have presented the sketch of simulations of relativis- 
tic plasma. We think that 10 5 — 10 7 processors are needed 
to perform the task. It requires much of intellectual efforts 
also, but it has to be done if we wish to figure out the 
microphysics of relativistic plasma. 



The issue still remains to be explained whether the 
plasma that allow for Cracow acceleration exists. The an- 
swer is simple, observations resolve the problem. Cracow 
acceleration predicts a lot of peculiar facts as superluminal 
shocks acceleration, acceleration at shocks with extremely 
large Lorentz factors (pulsars), the negative GRB lag, the 
extremely short acceleration time and others. 

Determining the acceleration time in pulsars would be 
the crucial proof. The plasma downstream of the shock in 
pulsars is at rest with respect to the interstellar medium 
(ISM) and the value of the homogeneous magnetic field can 
be detected there. The acceleration time in the mechanism 
of Cracow acceleration is equal to 1 or is shorter if measured 
in units o f the homogene ous magnetic field downstream of 
the shock (|Bednarzl . l2000t) . Mechanisms in which the change 
of the direction of the movement of a particle enabling it 
the return to the shock takes place downstream of the shock 
can not achieve such short acceleration time. 

3. The mechanism of Cracow acceleration 

In t his paragr a ph we follow simulations of lBednar3 ([2000) 
and iBednarzl (|2004|) . We will sometimes denote '(the 
plasma) upstream of the shock' as UP and '(the plasma) 
downstream of the shock' as DOWN. We have simulated 
relativistic shocks with Lorentz factors 7> 1 and the per- 
pendicular homogeneous magnetic field. Downstream of the 
shock the pure homogeneous magnetic field is present, and 
upstream of the shock there are additional magnetic field 
fluctuations. The magnetic field fluctuations pattern UP is 
characte rised by a value of the parameter Q that fulfils the 
formula (Bcd narzL l2004t) 

K±_ 1 

^ = 1 + Qp 2 ' 

where p is the particle momentum and k±, ku are perpen- 
dicular and parallel diffusion coefficients, respectively. 

The numerical calculations involve particles with mo- 
menta systematically increasing over several orders of mag- 
nitude. In order to avoid any energy dependent systematic 
effect we consider the situation with all spatial and time 
scales - defined by the diffusion coefficient, the mean time 
between scatterings and the shock velocity - proportional 
to the particle gyroradius r g , i.e. to its momentum. It means 
that the results are momentum independent and can be eas- 
ily scaled to any momentum. In particular we will be using 
the parameter Qq = Qp 2 that determines magnetic field 
fluctuations UP independently of the value of the momen- 
tum. 

We have applied the light velocity, c, equal to 1 and 
the time is measured in r g /c units. The plasma UP is hot 
and the compression is equal to 3U 2 , where U s is the shock 
velocity in the upstream plasma rest frame. The Lorentz 
factor of the plasma DOWN as seen UP, j u d, can be ap- 
proximated through the formula = O.7I7. 

Lengths of the magnetic field fluctuations required by 
Cracow acceleration must be much smaller than the mean 
time the particle spends UP between shock crossings, At u , 
where At u is measured in the plasma rest frame. It means 
that the time period in which the particle and a fluctuation 
interact is much smaller than At u . For better understand- 
ing we give At u in radians, for such an angle the particle 
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A 3 



Fig. 3. The angle <f> is increasing during the movement of 
the particle. The particle is entering UP with parameters 
from the set A3, and is leaving UP with parameters from 
the set A\ or A%. 

will turn in the wandering in the homogeneous field. The 
value is the same as in r g /c units. 

The mean change of the particle momentum direction 
between the moment when the particle is entering UP and 
the moment when it is leaving UP is equal to Aa u if mea- 
sured in the plasma rest frame and Acid if measured in the 
downstream plasma rest frame. Since upstream of the shock 
particles momenta are directed almost parallel to the shock 
velocity we assume that Aad ~ 2j u dAa u — lA2^Aa u . 

Let us consider the particle that crosses the shock from 
UP to DOWN (UP-DOWN) with parameters and r, 
where the particle phase £ (0, 2ir) is the angle between 
the direction DOWN-UP and the projection of the particle 
gyroradius onto a plane that is perpendicular to the homo- 
geneous magnetic field, r, see lBednarz fc Ostrowskil (|1999h . 
The sense of the angle is the same as the sense of rotation 
of the particle, and r, are measured in the downstream 
plasma rest frame. 

The particle will cross the shock again if there is a so- 
lution of the equation 

,t t t 

rsm(- + 0)sin- + - = 0, 

at positive time t. 

We will therefore distinguish three ranges of particle 
parameters (see Figure [3]), Ax = {0 : e (0i,02)} - the 
equation has not the positive solution and the particle does 
not return to the shock, A 2 — {(f) : € (02,03)} - the 
equation has the positive solution and the particle returns 
to the shock, A 3 = {0 : e (0 3 ,0i)} - the particle with 
parameters from this range crosses the shock DOWN-UP 
and can not cross the shock in the opposite direction, r sets 
the amount and determines values. 

Values of 0i and 03 are given by analytical formulae, 

01 = 2tt — arcsin(l/(3r)), 3 = 3ir/2 — arccos(l/(3r)), and 

02 must be calculated numerically. A particle with r < 1/3 
does not return to the shock. Chosen values of parameters 
are r = l,0i = 5.94, 2 = 1-96, 3 = 3.48; r = 0.5, 0i = 
5.55, 02 = 2.96, 3 = 3.87; r = 1/3, 0i = 2 = 3 = f. 

The change of the particle direction, Aa u , consists of 
two parts, the change in the process of wandering in the 
homogeneous magnetic field, Aw/,, and the change as the 
result of collisions with fluctuations, Awj, so we can write 
Aa u = AtJh + Aw/. We have neglected addition of vectors 
in this Section. 




ln(Q ) 



Fig. 4. The energy spectral index of accelerated particles 
versus Qo for shock Lorentz factors in the range 4-32768. 

A particle is entering UP having parameters (r, 0) from 
the set A3. If there are no fluctuations UP, the particle 
is returning DOWN having (r, 0) from the set A\. That 
is the case of Fermi acceleration, where the acceleration is 
forbidden in superluminal shocks and Aa u = AuJh- 

In Cracow acceleration Auif plays the key role. 
Fluctuations must be heavy so that Aujj is large enough. 
The angle change - Alu^ is directing a particle this way so 
that it gets parameters from the set A\ and Auj f is directing 
randomly, enabling the particle very much often changing 
from A3 to A2 through 03 . When Aw f dominates then half 
of the particles change to a value from the set A2 and 
half from the set A\. It implicates the short acceleration 
time equal to 1. 

The change of the particle direction through Aw/j is the 
linear process because At u is small for large 7, and through 
Aujf is the diffusion process, and that is expressed by the 
formulae (Awj) 2 — CfAt u and Aw/, = At u , where C/ is 
constant. It yields the relation Aa^ = \A2 r y(^jC '/ At u + 
At u ). If 7 grows then At u decreases. One can see now that 
the relative Awj contribution to Aad grows with 7, which 
means that the particle energy spectral index, /?, decreases. 
The index would increase but only if At u decreased more 
quickly than I/7 2 , but it does not occur as is shown below. 

Let us return to the simulations. We have performed 
simulations for 14 7 values, 7 = 2", where n — 2,3, ..,15 
and 500 Qo values from 5 • lO^ 1 to 5 • 10 7 in logarithmic 
space. For larger 7 we started the computations from larger 
Qo than 5 ■ 1Q _1 , but being small enough in order to get 
the satisfyi ng range of the asymptotic energy spectral index 
00 = 2.23 (|Bednarzl . l2004h . 

The obtained /3 — /3rj values as a function of ln(Qo) 
are shown in Figure H] Individual curves represent con- 
stant 7 values and one can see that they are similar. 
We have converted the curves according to the formula 
W = ln(Q ) + 2(ln(2.145)-ln(7)). After the transformation 
the curves have become identical, but the one for 7 = 4, and 
the curve for 7 = 8 has been standing out from this relation 
very slightly (see FigureJS]). The relation 7 2 /Qo = const for 
constant ft is exact and there must be a simple theoretical 
explanation for it. 

In order to get errors of fitting we have limited the 
range of the data to IF € (—4,2.5), see Figure [5] We 
have excluded points for which < /3q and curves for 
which 7 = 2 14 ,2 15 because their values do not cover 
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7 

6 
5 
,4 
' 3 
2 
1 

°-4 -3 -2-10 1 

ln(Q )+2(ln(2.145)-ln(Y)) 

Fig. 5. The same curves as in Figure HI but converted and 
cut for the purpose of fitting. The curve for 7 = 4 (dotted 
line) is standing out a bit from the fitting and for 7 = 
8 (dashed line) very little. Other curves fulfil the fitting 
closely and are represented by the curve for 7 = 256 (solid 
line). 





ln(Y"/Q ) 

Fig. 6. The values of At u j and Aa u j vs. ln(7 2 /Q ). There 
are two sets of 14 curves each for constant 7 = 2", where 
n = 2,3, ..,15, in the figure. Curves for 7 = 4 stand out 
a bit and are represented by dashed lines. Other curves 
are merging into one curve for Aa u and one for At u . One 
can see from the figure that for large 7 the value of Aa u 
decreases linearly with 7 and At u faster than that. 



the upper range of W . We have been fitting the relation 
P - p = Aexp(BW) and have gotten A = 1.007 ± 0.004, 
B = 0.722 ± 0.003. Finally, we have received the important 
relation 

(/3-/3o) 1 ' 39 = g (^) 2 , 7>8. (1) 

The numerical precision prevented us from carrying out 
calculations for 7 larger than 32768, but we can not see an 
obstacle that the above formula is not valid for all 7 > 8. 

We have repeated the above procedure for accurate val- 
ues of jud instead of 7 and have gotten the same results 
but the values of constants. Two first curves (7 € {4,8}) 
the same diverge from the relation as for 7. The separation 
of the curve for 7 = 4 depicts where Cracow acceleration is 
starting weakening, for 7 = 2 it turns off practically. 

Further, we have checked how At u and Aa u are chang- 
ing with 7. The choice of ln(7 2 /Qo) as independent variable 
enabled us to put curves in one line. In Figure [5] 14 curves 
for constant 7 = 2 n ,n = 2, ..,15 and in Figure [7] 14 curves 
for constant Qo € {1.4, 30, 297, .., 8.3 • 10 5 } are presented. 
The curves for 7 = 4 (dashed lines) and Qo — 1-4 are di- 
verging from general relations. 

From Figure [5] it is seen that Aa u j is constant for large 
7 and At u decreases with 7 faster than I/7. Moreover, we 
have put lines that are indicating values of P for fixed 7 2 /Qo 
and now one can estimate how reducing the contribution 
of Aujh to Actd is decreasing P to Pq. 

The aim of making Figure [7J was to show how At u was 
changing with 7. It decreases with 7 always slower than 
1 /7 2 and approaches the asymptotic value of 1 /7 2 for large 
7. It guarantees that Cracow acceleration does not fade 
with growing 7. For large 7, we have fitted the relation 

At u7 2 /V^o = 1.097 ±0.002. (2) 

Cracow acceleration could weaken with growing 7 be- 
cause of the increasing energy of reflected particles, but 
it is not. Let p s be the momentum of a seed particle. 
Reflection process increases particle energy ~ 7 2 times, so 



— 

"H 0,8 
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Q =1.4 ^ 








=30 J 


286<Q <S.3E5 












ln(r/Q ) 



Fig. 7. The value of At u <y 2 /y/Qo~ vs. ln(7 2 /Q ). There are 
14 curves for constant Qo in the figure. Curves for Qo = 1-4 
and Qo = 30 are indicated and other curves are merging. 
It is seen that At u j 2 / y/Qo is constant for large 7. 

that Qo ~ Q{Psl) 2 at the onset of Cracow acceleration pro- 
cess, so l/(Qp 2 ) ~ 7 2 /Qo = const if P = const. It means 
that the initial value of P depends on plasma conditions 
upstream of the shock only (Q and p s ). In real physical con- 
ditions this important result, similarly as previous results, 
can be altered by magnetic field fluctuations downstream 
of the shock. 

Cracow acceleration could weaken with growing 7 as a 
result of reducing the effectiveness of reflection process, but 
it is a quite different problem for which plasma conditions 
DOWN should be examined. 

4. Simulations 

Below, the light velocity is used as the velocity unit, c = 1. 
As the considered particles are ultrarelativistic ones, p = E, 
we often put the particle momentum for its energy. 

In the simulations we consider the mean magnetic field 
configuration perpendicular to the shock normal and mag- 
netic field fluctuations. We restrict our consideration to the 
test-particle approximation in which it is assumed that par- 
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tides are scattered by scattering centres in the fluid but 
have no effect either on the fluid velocity or on the density 
of scattering centres. Between two successive scatterings 
the particle is assumed to proceed along the undisturbed 
path in the mean field. We model particle trajectory per- 
turbations by introducing small-angle ra ndom momentum 
scattering along the mean field tra j ectory (lOstrowskil Il99lt 
iBednarz fc Ostrowski 119961 119981: [Bednard. l2004ft . 

The perturbed magnetic field represents the traditional 
picture based on the concept of magnetic scattering centres. 
It is simulated by the small amplitude particle momentum 
scattering within a cone with angular opening AD, less than 
the particle anisotropy ~ 1/7, where 7 is the Lorentz factor 
of the shock. The particle momentum scattering distribu- 
tion is uniform within the cone. For each Af2 the adequate 
mean time between scatterings was chosen to keep the same 
magnetic field perturbations pattern. 

The pattern is characterised by a value of the pa- 
rameter Q defined in Section [3] In this Section we use 
more convenient parameter to work with. We have defined 
C = — log 10 Q. If C determines magnetic field fluctuations 
upstream of an ultrarelativistic shock moving with a con- 
stant velocity, then with growing C particles are able to 
accelerate to higher energies in Cracow acceleration. 

Particle momenta are measured in the unit of any mo- 
mentum pq. We have applied seed particles momenta equal 
to 0.1 because of computational reasons, so that the value 
log 10 (p) = 1 means that the seed particle momentum was 
increased hundred times. 

The value of the mean magnetic field taken in the up- 
stream plasma rest frame, B u , is the unit of the magnetic 
field. The particle charge, eo, is the charge unit, so that the 
unit of time is po/ (B u eoc). 

Our simulations are intended to model decelerating ul- 
trarelativistic shocks, so that we have chosen the equation 
of motion of the shock in the form 



on, but for simplicity its speed is always equal to the speed 
of light. The rule yields the fields 



x s (t) 



-a In 



a + t 



where a > is a constant and x s , t are the distance and the 
time measured in the upstream plasma rest frame, respec- 
tively. 

There is a lack of knowledge of the plasma conditions 
downstream of the decelerating relativistic shock, so that 
we have defined them. We have divided the motion of the 
shock into discrete instants. We have applied a constant 
shock velocity at each instant. The downstream plasma ve- 
locity and magnetic field are derived from the consta nt ve- 
locity relativistic shocks (jBednarz fc Ostrowskil 119961 ) with 
the compression Rh = 2>U 2 for the hot plasma, where U s 
is the shock speed in the upstream plasma rest frame. The 
shock is perpendicular, so that the instant formulae are 



U(x,t) 



3t/ s 2 (t- At) - 1 
2U s (t-At) ' 



B(x,t) = -U a (t-At)* 



hU 2 {t-At)-l 
l-U?(t-At) ' 



(3) 



(4) 



where At is derived from At = x s (t — At) — x and x, t 
are the distance and the time respectively, x s is the shock 
position and all quantities apart from B(x,t) are measured 
in the upstream plasma rest frame. 

We were injecting mono-energetic seed particles up- 
stream of the shock with momenta distributed uniformly on 
the sphere. The seed particle density along the shock path 
was constant. Each particle trajectory was followed using 
numerical computations until the particle escaped through 
the free escape boundary placed far downstream from the 
shock or it reached the time larger than the assumed up- 
per limit. When the particle was wandering downstream of 
the shock its momentum was Lorentz transformed to the 
respective plasma rest frame at each small step and its tra- 
jectory was computed in the frame. 

The initial and final Lorentz factors of the shock were 
chosen. The factors and the equation of motion put the 
upper limit of time. We have divided the wandering time, 
measured in the upstream plasma rest frame, into ten equal 
periods. The successive time bins are labelled numbers from 
to 9. The Lorentz factor range is labelled jd = 7i~7/j 
where 7,, 7/ are the initial and final Lorentz factors, re- 
spectively. For example, 7^ = 40_5 means that we started 
collecting data at the Lorentz factor of the shock equal to 
40 and stopped at 5. 

The particle crossings of the shock were divided into 
stages as described in Figure U3 The stages are labelled 
'su', 'rd', 'ru', 'pi', 'pu', 'pd'. The spectrum of the particles 
crossing the shock from DOWN to UP for the first time is 
labelled 'rd', for the second 'pi', for the third and following 
times 'pd'. The total spectrum of the particles crossing the 
shock from DOWN to UP is the sum of 'rd', 'pi', 'pd' 
spectra and is labelled 'ad'. We introduced the token sp 
also, for instance sp='rd'. 

Each time a particle crossed the shock the respective 
contribution was added to the given momentum bin for the 
appropriate time bin and stage. The momentum and time 
are measured in the upstream plasma rest frame. 

We assume that particles emit synchrotron radiation. 
We neglect the self-absorptio n and take into ac count the 
total emission coefficient only (jPacholczvkL Il970h 



U = 



2>U 2 S - 1 
2U S 



B = -B U U S 



I9U? 



1 



1 — 17? 



S v = c 3 H sin 1? / N(E)F(x)dE, 
'0 



where U is the downstream plasma velocity close to the 
shock measured in the upstream plasma rest frame and B 
is the mean magnetic field close to the shock taken in the 
downstream plasma rest frame. 

We had to accept a rule in order to get U(x,t) and 
B(x, t) fields downstream of the shock. The rule is that the 
decelerating shock 'puts the instant fields back' and goes 



where C3 is a constant, H magnetic field intensity, $ the 
angle between the mean magnetic field and the direction 
toward the observer, N(E) - the particle energy distribu- 
tion, F(x) = xj K 5 / 3 (z)dz, K 5 / 3 (z) is the Bessel func- 
tion of the second kind, x — v/v c , v - photon frequency, 
v c = c\H siwdE 2 , c\ a constant and E is the particle en- 
ergy. 
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Crossing UP-DOWN 
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Stage 2 

DOWN 




UP 



pi 
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DOWN 




UP 



pd 



Crossing DOWN-UP 



Fig. 8. Stages of particle crossings of the shock, where 
DOWN is downstream and UP is upstream of the shock. 



We have calculated the spectral distribution of syn- 
chrotron radiation of accelerated particles using the formula 



N(p)F(-)dp, 



P 



where I(v) is intensity of radiation and N(p) the particle 
momentum distribution. Units of I(v) and v are unimpor- 
tant here and we applied a constant equal to 1. 

In order to represent each set of spectra concisely we 
have introduced the vector (7^, a, C, sp). For instance, 
(20_5, 10 2 , 3, pd) represents ten spectra of particles pro- 
duced by a relativistic shock slowing down from Lorentz 
factor 20 to 5 with the parameter of the equation of motion 
equal to 10 2 , the upstream magnetic field fluctuations pat- 
tern measured by C — 3 and particles crossing the shock 
from DOWN to UP for the third and following times. 

Fluctuations of the magnetic field downstream of the 
shock are unimportant now. They should be strong enough 
to allow effective reflections of seed particles. We have ap- 
plied there Cd — — log 10 50, where Cd is the value of C 
downstream of the shock. The parameter Cd is a variable 
and C a constant in Section [3 



5. Slowing down shocks 

The velocity and magnetic fields (Equations [3J SJ) down- 
stream of the shock change the energy of wandering parti- 
cles, so that they give some contribution to the total accel- 
eration. We did not analyse the amount of the contribution, 
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the particle will increase : 
~ 2 times of its energy : 
if .it returns UP. 



V = const 
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Fig. 9. In constant velocity ultrarelativistic shocks particles 
double their energy along every lap. 
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Fig. 10. In slowing down ultrarelativistic shocks particles 
spend more time upstream of the shock and receive larger 
momentum gains. 



but it is of little significance to the paper. As a matter of 
fact our main results do not depend on a kind of accelera- 
tion but it has to be fast enough. Moreover, the presented 
figures suggest that the contribution does not influence re- 
flection and Cracow acceleration processes significantly. We 
discuss the processes below, but one should remember that 
for strong deceleration they are affected by patterns of ve- 
locity and magnetic fields downstream of the shock. 

Cracow acceleration is very sensitive to the physical con- 
ditions upstream of the shock. Its action in shocks moving 
with a constant velocity is very well known. The particles 
double their energy along every lap there (Figure ^ and 
can not be accelerated if there is a lack of magnetic field 
fluctuations upstream of the shock. 

In slowing down shocks other factors give their contri- 
bution, but it is still the same mechanism since the acceler- 
ation depends mainly on the fluctuations upstream of the 
shock. 

Particles spend a bit more time (in units of the particle 
gyroradius) upstream of slowing down than upstream of 
constant velocity shocks. This effect has two consequences. 

First, particles UP have smaller anisotropy in slowing 
down shocks. That means that the particles receive larger 
energy gains along every lap (Figure [TU|) but their proba- 
bility of returning to the shock is smaller and in result they 
show steeper spectra and shifted to a higher energy. This 
effect is valid for particles that enter UP with momenta 
more parallel to the shock normal. Such particles increase 
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their momentum more and therefore the steeper and shifted 
spectra are seen at higher momenta. 

Second, particles are staying longer upstream of the 
shock in an absolute unit of time (Figure [10]). It means that 
particles that have higher energies return to the shock later. 
This entails that the delaying particles are stopping the 
high energy part of the spectrum from moving back to lower 
energies, what means that the part is moving back slower 
or standing or moving ahead. We call the phenomenon lin- 
gering. 

There is still an issue to be addressed. Some particles 
enter UP with the directions of the momentum close to the 
limit of the transition. Their interaction with magnetic field 
fluctuations UP is weak and their increase of the momen- 
tum is small. If the shock has a constant velocity then a 
fraction of the particles can not reach the shock again, but 
when the shock is slowing down the particles get an addi- 
tional time and change the trajectories in such a way that 
some particles from the fraction can reach the shock again. 
This leads to flatter spectra at lower energies. 

In result, Cracow acceleration in slowing down shocks 
produces flatter spectra at lower energies and steeper at 
high energies. 

The stronger deceleration of the shock the shorter dura- 
tion of the episode if "fi and 7/ are fixed. It means that the 
total time of the episode could be very small what yields 
two effects. 

The first effect refers to reflection process. Reflecting 
particles gain higher energies if they are spending more 
time downstream of the shock. At first to the slowing down 
shock (from DOWN) arrive particles which increase the mo- 
menta less and then with better energy gains. These all 
particles have entered DOWN at the same Lorentz factor 
of the shock. This gives smaller range of energies (across 
all time bins) of particles reflecting from strong than from 
weak decelerating shocks. 

The second effect refers to Cracow acceleration. 
Particles accelerating at strong decelerating shocks receive 
smaller energy gains because they have not enough time to 
receive larger. 
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6. Results 

From graphs of particle spectra one can get the energy spec- 
tral index 

, , _ A[log 10 (dJV)/dlog 10 (p)] 
A[log 10 (p)] 

In each of photon spectra we distinguish the frequency 
Vmax where the maximum intensity I{v max ) is reached. In 
principle, we will limit ourselves to frequencies larger than 

The first lap in the process of Cracow acceleration is rep- 
resented by 'pi' spectra. They are similar to corresponding 
'pd' spectra. We have not presented a 'pi' spectrum, but 
the observation is visible as one compares Tables \T\ [5] and 
|U Our first result is that reflection process does not influ- 
ence Cracow acceleration from the first lap, and therefore 
we can treat these processes as separate. 

In Figures QTJ H2J H3] and [14] are presented spectra that 
are produced by practically constant velocity shocks. The 
flat part of 'pd' spectra (Figure [2]) has got j3 in the range 
2.5-2.7 what is larger than the limiting value of f3 produced 
by Cracow acceleration in constant velocity shocks. 

In next three Figures (ITBl HM ITT1) the particle spectra 
produced by a strong decelerating shock are shown. If they 
are compared with the three preceding figures then one can 
see to what extent some remarks presented in the previous 
Section are correct. 

Our primary task is to present how lags are formed in 
photon spectra of slowing down ultrarelativistic shocks. We 
have examined the spectra and collected results in tables. 
The description of each table contains the values of 7^, sp 
and two first columns the values of a, C . 
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Table 7. 7^ - 20_5, spectrum - ad 
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A lag will be visible as an intersection of two consecu- 
tive photon spectra (produced in two consecutive periods of 
time). Let as denote the frequency at which the intersection 
occurs as Vi , the number of the first spectrum as n and the 
second as n + 1 . If photon intensity of the n + 1 spectrum is 
higher than the intensity of the n spectrum just before i>i 
(at lower frequency), then the positive lag appears. If the 
intensity of the n + 1 spectrum is lower than the intensity 
of the previous spectrum just before z/j, then the negative 
lag appears. 

The first occurrence of the positive lag in the set of 
photon spectra for consecutive periods of time was put into 
columns third and fourth, where n p is the number of the 
previous spectrum and v p is the frequency at which the 
intersection occurs. The first occurrence of the negative lag 
was collected in the fifth and sixth columns, where n n is the 
number of the previous spectrum and v n is the frequency 
at which the intersection occurs. 

The seventh column is marked by 'lag' and contains the 
information how lags in the photon spectra are changing. 
The photon spectra that are only increasing (decreasing) 
their intensities at all frequencies are denoted as R (D). 
If the intensities at first are increasing and next they are 
decreasing at all frequencies, then they are denoted as RD. 

If R appears at some values of a and C, then for all 
smaller values of a and for the same value of a and smaller 
values of C R appears also. In that case the subsequent 
items of R are not included in the tables. If D appears at 
some values of a and C, then for all larger values of a and for 
the same value of a and larger values of C D appears also. 
In that case the subsequent items of D are not included in 
the tables. 

If n spectrum and n + 1 spectrum build a positive lag at 
frequency v v and then n + 1 spectrum and n + 2 spectrum 
build a positive lag at frequency v p 2 < v v and so on (or 
not), then the spectra are denoted as PJ,. If only one pos- 
itive lag appears, then the spectra are denoted as P. Here 
it is necessary to add that we did not take into account 
positive lags that appeared at frequencies lower than v max . 
In fact positive lags are always Pi and are moving very fast 
towards low frequencies passing through v max . 

If n spectrum and n + 1 spectrum build a negative lag 
at frequency v n and then n + 1 spectrum and n + 2 spec- 
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trum build a negative lag at frequency v n 2 < v n and so on 
(or not), then the spectra are denoted as Ni. If only one 
negative lag appears, then the spectra are denoted as N. 
The negative lag is moving slowly towards low frequencies 
and stops at high frequencies far before v max . 

If n spectrum and n + 1 spectrum build a negative lag at 
frequency v n and then n + 1 spectrum and n + 2 spectrum 
build a negative lag at frequency v n i > v n and so on (or 
not), then the spectra are denoted as Nt- If only one nega- 
tive lag appears, then the spectra are denoted as N. These 
types of spectra are always Nt- The case of N arises due to 
the limitation of data and large errors at high frequencies. 

An exception exists for small initial Lorentz factors of 
the shock and small a (Figures [18 1 [19 ] I2"0 J) . In Table[3]for the 
case of (10_5, 60, any, ad) the spectra number 8 and 9 build 
the negative lag N (FigureHH]). First, we do not expect such 
small value of a for small initial Lorentz factors. Second, a 
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Fig. 11. Photon spectra of (40_5, 10 4 , 3, ad). 




Fig. 12. Particle spectra of (20.5, 10 5 , 3, ad). The flat part 
of the spectra has got j3 in the range (2.4, 2.5). 
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Fig. 13. Particle spectra of (20_5, 10 5 , 3, rd). 



must grow from a moment. We expect that in real physical 
conditions such shocks form the spectra of the type Nf. 

The structure PJ, often appears together with the struc- 
ture NX, in this case we will consider the structure P^NJ,. 
The untypical structure PNf in (40_5, 10 3 , 3, ad) arose 
from the small time resolution. An example of the struc- 
ture PJ,NJ, without details is presented in Figure [21] and 
with details in Figure |2"2"1 

Data presented in tables with 'rd' (Tables [4l [5j and 
'pd' (Tables Q] [5]) spectra give us the very important re- 
sult that individual processes give their own contributions 




i°g, (p) 



Fig. 14. Particle spectra of (20_5, 10 5 , 3, pd). The flat part 
of the spectra has got /3 in the range (2.5, 2.7). 




log,„(p) 

Fig. 15. Particle spectra of (4CL5, 60, 3, ad). Maximum 
momenta are small because particles have not enough time 
to get larger. The value of j3 of the flat part of the spectrum 
at the time 2 is in the range 1.2 — 2.2, at the time 4 is equal 
to 1.0 and at the time 7 is equal to 0.7. 




Fig. 16. Particle spectra of (40_5, 60, 3, rd). The smaller 
range of maximum energies of particles reflecting from the 
strong decelerating shock is visible. 



to the general spectrum. Reflection process generates the 
structure P| (the positive lag) and Cracow acceleration the 
structure Nf (the negative lag). That is the only result we 
have obtained from the tables. Henceforth, the tables with 
'ad' (Tables El El [U El US spectra will be discussed. 
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Fig. 17. Particle spectra of (40_5, 60, 3, pd). Cracow ac- 
celeration process in slowing down shocks produces natter 
spectra at lower energies and steeper at high energies. 




i°g 10 (p) 

Fig. 20. Particle spectra of (10.5, 60, 2, rd). The high- 
energy tail moves back too slowly to be able to produce 
a positive lag. 
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Fig. 18. Photon spectra of (10_5, 60, 4, ad). In real physical 
conditions the shock should produce the lag of the type Nf. 
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Fig. 19. Particle spectra of (10_5, 60, 2, ad). The value of 
f3 of the flat part of the spectrum at the time 3 is equal to 
1.7, at the time 4 is equal to 1.4 and at the time 7 is equal 
to 0.9. 
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Fig. 21. Photon spectra of (20_5, 60, 4, ad). An example of 
the structure P|N| without details. 
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Fig. 22. Photon spectra of (40_5, 10 2 , 3, ad). An example 
of the structure P|N| with details. 



When one is examining the tables, one can see that the 
structure PJ-NJ. occurs at large initial Lorentz factors of the 
shock and strong decelerations. If % is smaller and/or a is 
larger then the spectra exhibit the structure Nj. 

To obtain a better resolution of (40_5, 10 3 , any, ad) and 
(20_5, 10 3 , any, ad) spectra we have performed the same 
simulations as previously but 7/ was equal to 10 instead of 



5. The additional simulations were following (40_10, 10 3 , 2, 
any), (40.10, 10 3 , 4, any), (20.10, 10 3 , 2, any), (20.10, 10 3 , 
4, any) and the outcome was put into Tables M and ITTfl 

The spectra of (40.10, 10 3 , 2, ad) show the structure 
P4- at low frequencies and the unimportant structure NJ, 
at higher frequencies and later times what is not visible 
for = 40.5. The spectra represent almost pure PJ, that 
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Fig. 23. Photon spectra of (40.10, 10 3 , 2, ad). The spectra 
are practically purely P^. 
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Fig. 24. Photon spectra of (20.10, 10 3 , 2, ad). We have 
labelled the spectra NNf since the point N does not move 
at first. 

is not suppressed by Cracow acceleration (Figure |2"5|) , we 
neglected results if v was high (here P for v — 4.75). 

The case of (40.5, 10 2 , 4, ad) is in fact P|NI, but N| is 
outside the range of detectable results (Table [5]). 

We have analysed the particle and photon spectra and 
reached the following conclusions 

a) the structure P\, is formed when the sudden fall in 
the spectra of reflected particles 'rd' is not suppressed by 
Cracow acceleration, 

b ) the structure P\.N\. is formed when the growth of the 
high-energy tail of Cracow acceleration spectra 'pd' is ac- 
companied by the sudden fall in the spectra of reflected par- 
ticles 'rd', 

c) the structure N\ is formed due to lingering of the 
high-energy tail of Cracow acceleration spectra 'pd' in the 
presence of the weak fall in the spectra of reflected particles 
'rd'. 

The fall in the spectra of reflected particles is best visible 
when one is looking at the high-energy tail of the spectra. 
Before doing so it is convenient to raise the particle spectra 
of initial times to the density level of the remaining spectra. 
Then, the high-energy tail will almost always move towards 
low energies. 

We have done rough estimation of the fall and intro- 
duced the parameter: —Ap/At, where Ap is the difference 
of the momenta measured at a constant particle density 
between the two particle spectra that differ by the time 
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Fig. 25. Particle spectra of (20.5, 10 3 , 1, ad). The lingering 
of the high-energy tail that arises from Cracow acceleration 
operating at the slowing down shock is visible. The value of 
fj of the flat part of the spectrum at the time is equal to 
1.2 and then there are two flat parts, at the time 2 /3 = 0.7 
and P = 2.1, and at the time 5 j3 = 0.9 and j3 = 1.8. 
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Fig. 26. Particle spectra of (20.5, 10 3 , 1, rd). 
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Fig. 27. Photon spectra of (20.5, 10 3 , 2, ad). An example 
of the structure Nf. 



At. The value of the parameter was measured at the high- 
energy tail. The results are collected in Table [TT] In the 
table ni, n2 are numbers of spectra for times t\, t2, where 
momenta p\ , p2 are measured at a selected density, respec- 
tively. 
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Table 11. Power of falling of 'rd' spectra 
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In spectra of (20.10, 10 3 , 2, rd) we have observed a weak 
structure PI at frequencies lower than v max . We have es- 
timated that the case is close to the boundary between 
appearing and not-appearing of Pi and the boundary is at 
— Apj At ~ 8. One should remember that both the param- 
eter and the measurement are rough. 

We denoted the spectra of (20.10, 10 3 , 2, ad) as NNf 
since N keeps still at first (Figure |24|) . 

The lingering of the high-energy tail in Cracow accelera- 
tion spectra is presented in Figure [231 and related 'rd' spec- 
tra in Figure[5U The tail moves ahead if C is larger. Photon 
spectra with the structure Nf are presented in Figure |2"T1 
The phenomenon of the production of lags is simple. 
When one of the two acceleration processes dominates the 
other then appears Pi or Nf. When they compete then Pi 
is pulling the progenitor Nf and changes it into N.J,. This 
interaction is mutual and the progenitor Nf is damping 
down Pj. what makes it less distinctive. 

Negative lags appear at frequencies higher than v max 
only and the difference between the frequencies and v max 
is not small then. Positive lags are found at frequencies both 
lower and higher than v max but mainly around v max . 

We normally expect that shocks with the small value of 
7i have large a and vice versa. If we adopt the above as- 
sumption then, according to the received results, there is a 
certain value of the initial Lorcntz factor 7$ below which the 
positive lag is not visible and with growing 7^ positive lags 
become more and more harsh. If we assume additionally 
that the GRB luminosity is growing together with increas- 
ing 7j, then together with the increase in the distance to 
GRBs we will be observing the larger contribution of GRBs 
with stronger positive lags. 

Our simulations show that the low energy electrons gain 
small f3 if the shock decelerates quickly and if 7 decreases 
then j3 decreases too (see Figures Ql}] and [TH]) . It is con- 
sistent w ith observations of the low energy electrons in 
lerions dWeiler fc Panagial fl978l 1.2 < /? < 1.6), GRBs 



P 
d 



Band et al. 19931 6 ^1.1 typically), hotspots in radio 



galaxies (IStawarz et al.ll2007l . /3 < 2) and luminous blazar 
sources (|Sikora et al.l [20091 /3 ~ 1.6 typically). Moreover, 
the simulations show that the high energy electrons gain j3 
larger than 2.23 if the shock decelerates slowly (see Figure 
[T2l 2.4 < j3 < 2.5) and it is consistent with observations of 
the sources with relativistic shocks where f3 ~ 2.4 is a typ- 
ical value. Medium deceleration is presented in Figure l25l 
there are two flat parts in spectra. Our intuitive estimation 
is that the deceleration of a typical GRB is changing from 
a = 50 to a = 9 • 10 3 during the phenomenon. 

7. Fermi acceleration 

In order to check whether Fermi acceleration is able to pro- 
duce a GRB lag we have carried out additional Simula- 
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tions. We have switched Cracow acceleration off through 
the application of C = —2 upstream of the shock for all 
simulations and we have switched Fermi acceleration on by 
choosing large values of Cd downstream of the shock. In 
this section we consider various values of Cd and we use 
the adapted vector (7^, a, Cd, sp) instead of the old one. 

At fixed values a = 10 2 and Cd = 5 we carried out sim- 
ulations for 7d G {10.5,20.10,40.10}, and at fixed a = 10 3 
and two values of 7^ G {10.5,20.10} we performed simu- 
lations for Cd G {1, 3, 4, 5, 6, 7}. It turned out that Cd = 1 
was a saturated value for the second set and increasing it 
did not change the results. The results have been collected 
in Tables IMil 

We have checked spectra 'rd' and 'pd' and it turns out 
that Pi is a pure outcome of reflection process because 
spectra 'pd' do not show a lag in this case (Figure |2"5]) . The 
negative lag arises in full measure from Fermi acceleration 
what is visible in spectra 'pd'. In this case spectra 'rd' are 
showing Pi which is suppressed by Nt from spectra 'pd' 
and it gives spectra 'ad' with Nt (see Figure |2"9")) . 

Findings are simple. Fermi acceleration generates nega- 
tive lags only in shocks that have small Lorentz factors and 
are slowing down slowly. In other conditions positive lags 
produced by reflection process are visible. 

The first question is which way the negative lag is pro- 
duced. We have not prepared an expert opinion on this 
issue, but it has emerged a simple interpretation after 
analysing particle spectra 'pd' (Figures l30l l31"j) . Spectra in 
Figure [301 are similar to the corresponding spectra with ac- 
tive Cracow acceleration. The spectra are falling at medium 
particle energies and moving forward at high-energy tail. 
In the process of Fermi acceleration particles are wander- 
ing UP over a distance shorter than in Cracow accelera- 
tion, but for small Lorentz factors of the shock it is enough 
to produce lingering of the high-energy tail in the particle 
spectra. 

The second question is why there is no negative lag for 
(10.5, 100, 5, ad). It results probably from the fact that 
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Fig. 28. Fermi acceleration. Photon spectra of (4CL10, 10 2 , 
5, ad). Fermi acceleration does not influence Pi produced 
by reflection process. 
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Fig. 31. Fermi acceleration. Particle spectra of (2CL10, 10 3 , 
3, pd). The corresponding photon spectra do not show nei- 
ther negative nor positive lags. 
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Fig. 29. Fermi acceleration. Photon spectra of (10_5, 10 3 , 4, 
ad). Fermi acceleration generates Nf and reflection process 
is too poor to suppress this. 
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Fig. 30. Fermi acceleration. Particle spectra of (f 0_5, 10 3 , 
1, pd). The spectra are similar to the corresponding spectra 
with active Cracow acceleration producing Nj. 

Fermi acceleration has longer acceleration time and pro- 
duces steeper spectra than Cracow acceleration. It explains 
the lack of negative lags in photon spectra of rapidly slow- 
ing down shocks. 

The third question is why one can not see negative lags 
for larger Lorentz factors of the shock. The answer is that 
the wandering time UP decreases more quickly with the 



growing Lorentz factor of the shock than the time in Cracow 
acceleration or Fermi acceleration fades out. 

On the basis of the received results we put forward the 
hypothesis that the negative lag can be produced by any ac- 
celeration mechanism in which particles after reflection pro- 
cess are crossing the shock repeatedly from UP to DOWN 
and around DOWN to UP. 

We have not examined Fermi acceleration in detail and 
we think that the Weibel instability and the fact that parti- 
cles are accelerating in PWNe make this mechanism impos- 
sible. Our numerical code is still relatively straightforward 
and therefore it can not be ruled out on the basis of the 
current results that the negative lag arises due to Fermi ac- 
celeration. More detailed simulations and observations are 
needed in order to decide conclusively on the matter. 

8. Two acceleration processes 

Observations of GRBs with the Fermi Gamma Ray Space 
Telescope indicate that photons in the GeV range (the 
LAT) and photons in the MeV range (the GBM) arrive 
at different times. The GeV photons arrive so metimes later 
and sometimes before t han th e MeV photons (jAbdo et all 
120091 : iGhirlanda et all l2010t ). It seems that the observed 
spectra in the two ranges originate from different regions. 
However, it is the very same phenomenon that occurs in 
one region and we will describe it below. 

In order to understand this phenomenon one should re- 
alise that Cracow acceleration does not have a general re- 
striction to the maximum energy achieved by particles. The 
restrictions result from specific physical conditions only. 
Values of C and 7 are two basic parameters that are lim- 
iting the energy of accelerated particles. Increasing any of 
them increases this limit. 

The acceleration time is another important quantity and 
although it is constant downstream of the shock, it is 7 
times smaller in the upstream plasma rest frame and grows 
together with reducing 7. 

The Fermi Telescope observes bursts which have suffi- 
ciently large values of C and 7 to allow for Cracow accelera- 
tion to take place at GeV energies detected on Earth. When 
the Fermi Telescope is detecting the onset of the GeV emis- 
sion then the accelerated particles achieve energies enabling 
them to produce synchrotron emission in this range. 
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The two acceleration processes behave differently during 
the phenomenon. Reflection process is weakening almost all 
the time except for the short initial time when the number 
of particles entering the process is larger than the number 
of particles coming out. The weakening manifests itself by 
the moving of the maximum of the energy distribution of 
accelerated particles toward lower energies. We will call the 
maximum the reflection surge. 

Cracow acceleration accelerates particles all the time 
but unevenly. The acceleration time decreases and this in- 
fluences each particle the same. Reflection process provides 
particles with energy that is more and more lower, there- 
fore Cracow acceleration accelerates more and more smaller 
number of particles at a fixed particle energy. That does not 
concern particles with the highest energies. 

We expect that the real equation of motion of the shock 
is much more dynamic than that in our simulations. The 
value of a increases with decreasing 7. The shock with a 
small value of a has more particles UP than the shock with 
a large value of a. During the rapid phase the shock pro- 
duces a large amount of high-energy particles UP. The par- 
ticles make up an additional fraction of non-thermal par- 
ticles UP. It is convenient to regard the fraction as seed 
particles produced in the rapid phase. We will call the frac- 
tion the rise-fraction. 

When the shock slows down and a increases then the 
Cracow acceleration takes up the rise-fraction. This pro- 
duces a weakening wave moving toward higher energies in 
the energy distribution of accelerated particles. We will call 
the wave the Cracow surge. 

The question is why we can not see the Cracow surge 
in the present simulations. It results from the simplicity of 
the equation of motion which does not allow for a large 
dynamics. The simulations show the lags because they are 
strong local phenomena. The Cracow surge is a global phe- 
nomenon, but there are some signs of it in the present cal- 
culations. 

One can see bumps at high energies in 'ad' particle spec- 
tra where the Cracow surge should occur. At later times 
the bumps move toward low energies, but it is also caused 
by the low dynamic range. The corresponding 'ad' photon 
spectra show an extended plateau instead of the maximum, 
but better figures have been not presented in the paper. In 
real conditions the rise-fraction is produced quickly (small 
a) and then Cracow acceleration processes it slowly (large 
a). 

A GRB starts at the energy of accelerated particles that 
is equal to the energy of reflected particles at the onset of 
the phenomenon. A short time after the onset two bulges 
appear in the energy spectrum, the reflection surge moves 
toward low energies and the Cracow surge moves toward 
high energies. The transition of a surge is seen as a signal. 
Each of two detectors of photons working in different spec- 
tral bands can detect one signal only. The sequence of the 
detections depends on the GRB. 

There is the option that an instrument could detect the 
reflection surge at the very onset of a GRB. Shortly af- 
ter the detection the same instrument detects the Cracow 
surge. The light curve produced by the instrument will ex- 
hibit the structure with two pulses (the double pulse). Such 
structures are present in BATSE data. 

We describe the Cracow surge in more detail below. We 
will u se the observation of GRB 090510 ()Ghirlanda et all 
l2010f ) as the example. 



GRB 090510 was observed at MeV energies at first. This 
pulse could originate both from the reflection surge as well 
as from the Cracow surge. The pulse at GeV energies orig- 
inates from the Cracow surge. The period of time between 
the pulses is considerable, therefore we can assume that 
during the GeV pulse the value of a is large and constant. 
The rise-fraction is present upstream of the shock during 
the pulse. 

The rise- fraction consists of particles that are more close 
to the shock and have lower energies and larger anisotropy, 
and farther with higher energies and smaller anisotropy. 
The anisotropy does not play a significant role in the case. 
The acceleration time changes very slowly for large a and 
we can treat it as constant. 

At first, the shock is getting the particles of the rise- 
fraction with lower energies and accelerates them. When 
the shock reaches the particles with higher energies, then 
the lower-energy particles were in time to accelerate to en- 
ergies comparable to the higher-energy particles. As a result 
instead of the value of j3 typical for given a, e.g. j3 = 2.6, a 
flatter particle spectrum is arising, e.g. /? = 1.4, which will 
approach /3 = 2.6 going to higher energies, but will enter 
the GeV photon spectrum range with /3 — 1.6 for exam- 
ple. Next, the normal seed particles, that produce /3 = 2.6, 
are falling within the GeV range and the photon spectrum 
starts falling. 

The above description is compatible with the observa- 
tion of GRB 090510. At the high-energy part of the photon 
spectrum the photon index is equal to j3/2 + 1. In the GeV 
range the photon spectrum at first grows with the photon 
index equal to 1.8 (j3 = 1.6), and then it is falling with the 
photon index equal to 2.3 (j3 — 2.6). 



9. UHECRs acceleration 

We claim that the mechanism of Cracow acceleration pro- 
duces the ultra-high energy cosmic rays (UHECRs). The 
restrictive requirement that rules out UHECRs acceleration 
in most astrophysical objects is that the particle gyroradius 
must be much smaller than the system size. The restriction 
does not exist if the acceleration takes place at standing 
or moving with a mildly relativistic velocity shocks (with 
respect to ISM) and UP has the Lorentz factor larger than 
~ 1000 (with respect to ISM) and UHECRs come out from 
UP. Two kinds of sources fulfil the requirements, PWNe 
and relativistic flows in AGN jets. The PWN sources over- 
come the restriction of the Greisen-Zatsepin-Kuzmin cutoff 
also. 

Below, we explain how Cracow acceleration is able to 
accelerate UHECRs. All quantities are given in the plasma 
rest frame. Let E u be the particle energy UP, Ed the parti- 
cle energy DOWN, B u the homogeneous magnetic field UP 
and Bd the homogeneous magnetic field DOWN. When the 
particle crosses the shock from UP to DOWN then its en- 
ergy changes according to the formula Ed = E u /-f u d- When 
the particle crosses the shock from DOWN to UP then its 
energy changes according to the formula E u = 2"/ u dEd- This 
is a result of the particle anisotropy upstream of the shock. 
The relation between magnetic fields is Bd = A-f u dB u . We 
assume that DOWN is at rest or moves with the velocity 
< 0.95c with respect to ISM. 

Synchrotron losses are ~ E 2 B 2 . It gives that the syn- 
chrotron losses of a particle that crossed the shock from 
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DOWN to UP are ~ 4 times slower than the losses be- 
fore the crossing. The particle spends much shorter time 
UP than DOWN so that the synchrotron losses UP can 
be neglected. Synchrotron losses of protons and nuclei are 
neglected DOWN and their achieved maximum energies 
DOWN depend on the size of the system only. 

The accelerated particles do not escape from DOWN if 
the perpendicular size of system is larger or equal to the 
particle gyroradius r g . It yields the formula for the maxi- 
mum energy DOWN, Ed — SOOddBd [eV], where dd [cm] 
is the size of DOWN perpendicular to the shock normal, 
Bd [Gs] . Protons and nuclei are able to reach energies equal 
to ~ 10 14 — 10 16 eU downstream of the shock and it is 
enough to produce UHECRs if the Lorentz factor of the 
shock 7 > 10 3 and the particles get out into ISM through 
UP. 

We assume that a particle crossed the shock from 
DOWN to UP and its momentum was parallel to the shock 
normal at the beginning of the wandering UP. After a while 
the direction of the particle momentum changes and a is 
the angle between the momentum and the shock normal. 
The angle a is small therefore the distance the particle cov- 
ered in the direction perpendicular to the shock normal is 
d u ~ r g a 2 /2, where r g is the particle gyroradius UP. The 
value of d u is the same both in the upstream and in the 
downstream plasma rest frame but the particle gyroradius 
grows ~ &lud ti mes when the particle crosses the shock 
from DOWN to UP. The perpendicular (to the shock nor- 
mal) length of the particle path DOWN is dd — r g and 
therefore d u /d d ~ 47 2 rf a 2 . 

We will estimate the value of a in two ways. First, for j3 
slightly larger then /3 , we apply a ~ At u . We eliminate Q 
from Equations ffl [2] and get At u ~ 0.51(/3 - 2.23) 695 /7- 
If we put the values /3 = 2.3, 2.4, 2.6 into the equation we 
get d u /dd ~ 0.013,0.045,0.13. 

The electromagnetic forces that disturb the particle 
movement in the homogenous magnetic field twist its path 
with different radii of curvature. If we exclude the contri- 
bution of the homogenous magnetic field then we get the 
change of the direction of momentum and displacement of 
the particle caused by a scattering and, in consequence, 
the scattering radius. We will use the parameter r\ that is 
the ratio of the particle gyroradius to the mean scattering 
radius. 

In our simulations rj is equal to infinity. In real plasma 
rj has a finite value. If 77 takes a finite value then the first 
estimation of d u /dd becomes larger for half of the particles. 
Because d u /dd ~ a 2 , the first estimation gives that the 
deceleration of the shock must increase a a few times to 
allow accelerated particles come out from UP. 

The second estimation is in the case of the neglected 
homogeneous magnetic field. The angle a results from scat- 
terings only and = /3q. We apply a ~ Aa„, where 
Att„ — 1-3/7 ( see Figure HJ, but we divide the result by 77 
because the particle turns in electromagnetic disturbances. 
It yields the second estimation d u /dd — 3. 4/77. We suggest 
that loop plasma has the value of rj in the range 10 — 100. 
It gives that if a increases ~ 2 times up to a few times 
because of the deceleration of the shock then the particles 
come out from UP. 

In the model of loop plasma the homogeneous magnetic 
field and electromagnetic disturbances are generated by 
electric current loops. The loops could change if 7 changes. 
We think that B u could decrease with 7. It would explain, 



at least partly, dissipation of the Poynting flux of the pulsar 
wind before it reaches the inner edge of the Crab Nebula. 

A few times increase in the value of a as the result of 
the shock deceleration is something expected. The slowing 
down velocity of the shock enables particles to spend longer 
time UP and it can not be neglected. If B u decreases then 
r g increases and, as a result, d u increases and the particle 
spends longer time UP. It could even produce feedback that 
prevents the particle from returning to the shock. However, 
the reduction of the shock velocity can not be too large 
because the energy Ed increases ~ Zjudlis times, where 7i S 
is the Lorentz factor of the plasma UP as seen in ISM at 
the moment when the particle comes out from UP. 

The presented calculations are promising but further 
numerical simulations are needed to show that Cracow ac- 
celeration accelerates UHECRs. The first kind of simula- 
tions will be similar to presented in this paper but the nu- 
merical code will include an extended equation of motion 
and variable B u . The second kind of simulations should 
search for values of 77. One must choose parameters of loop 
plasma and follows particle trajectories in its electromag- 
netic field. The third kind of simulations will be similar to 
the first kind but UP will be replaced with loop plasma. 

If Cracow acceleration accelerates UHERCRs then as- 
tronomical observations will show that UHERCRs arrive 
from systems where the shock move away from an observer 
and is slowing down from 7 > 10 3 to 7 > 10 2 . 

Fermi acceleration can not accelerate UHECRs because 
its change of the direction of the particle momentum that 
enables the return to the shock takes place DOWN. 

10. Seed particles 

In our model presented in iBednarzl (|2004fl we have pro- 
posed some acceleration mechanisms to be responsible for 
production of seed particles and that their radiation is de- 
tected as the precursor of GRBs. Actually, the mechanisms 
arc ineffective in prc-accclcration of protons and nuclei and 
the analysis of spectra of GRB precursors is suggesting 
that they are produced in the same mecha nism as GRBs. 
In our opinion, the mechanism described by iHoshino et al.l 
(|1992l ) is the only (among the proposed) effective mecha- 
nism that is able to produce seed leptons. This is the pro- 
cess of the downstream plasma acceleration of positrons to 
non-thermal distributions. Cracow acceleration needs seed 
particles upstream of the shock to accelerate particles. We 
propose two ways in which they could be injected into UP. 

The first way occurs in the relativistic jets that are fired 
from a massive object close to the neighbourhood. The jet 
production mechanism and indeed the jet composition on 
very small scales are not known at present. We propose 
that seed particles are present in jets from the beginning 
of producing them and are there all the time. We think 
that loop plasma generates seed particles both leptons and 
nuclei. 

The second way is at work in a pulsar wind when it is 
streaming into the ambient medium and creating a stand- 
ing shock. Seed particles can, similarly as in jets, be present 
from the moment of generating the wind, however it con- 
cerns leptons only. We propose that seed particles, particu- 
larly protons and nuclei, are entering the pulsar wind from 
its side boundary. These particles are coming out from the 
ambient medium being pushed by the gravitational force 
of the pulsar. Their Lorentz factor in the wind plasma rest 
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frame is equal to ~ y w , where j w is the Lorentz factor of 
the pulsar wind, so that they become seed particles auto- 
matically. 

11. Summary and Discussion 

Up to now, simulations of particle acceleration at relativis- 
tic shocks concerned the constant velocity shocks only. We 
are the first who have simulated particle acceleration at 
slowing down relativistic shocks and that is why we have 
solved the mysteries of photon and particle spectra of the 
astronomical objects that radiate the synchrotron radiation 
of particles accelerated at relativistic shocks. No one was 
trying to perform such simulations before because it was a 
very difficult task. It turns out that GRB lags and the hard 
low energy spectral indices arise due to the deceleration 
(Section E]). 

We have presented the mechanism of Cracow accelera- 
tion (Section[3]). The mechanism increases its strength when 
the Lorentz factor of the shock grows and initial j3 does not 
depend on the Lorentz factor but on the conditions of the 
upstream plasma. 

We have shown that the whole acceleration consists of 
two processes, reflection process and Cracow acceleration. 
We have assumed that the real deceleration is more dy- 
namic than that in our simulations and the rise-fraction 
of particles upstream of the shock if formed at the strong 
deceleration stage. Acceptance of the assumption gives ex- 
planation to the question why two transitions are observed 
at one GRB and why the observation at high energy can fol- 
low the observation at low energy and vice versa (Section 
[8]). It also explains why the photon spectrum is hard at 
the rise stage of the Cracow surge. Sometimes, two surges 
produce the double pulse, the first maximum in the pulse 
originates from the reflection surge and the second from the 
Cracow surge. 

Cracow acceleration requires strong disturbances of the 
movement of particles upstream of the shock in order to 
accelerate the particles. We have shown that current PIC 
simulations of relativistic plasma are useless. We have pre- 
sented the model of the microphysics of relativistic flows 
that would allow for acceleration in the mechanism of 
Cracow acceleration (Section^). We have sketched the PIC 
simulations that are needed to solve the problem of the mi- 
crophysics of relativistic plasma. 

The current PIC simulations suggest that seed particles 
could be produced in the upstream plasma. We agree with 
the idea, but only proper PIC simulations can prove this. 
We proposed that particles of the ambient medium are in- 
jected into pulsar wind of PWNe and become seed particles 
(Section [Til. 

We claim that Cracow acceleration produces UHECRs 
which come out of UP (Section [9]). Relativistic shocks are 
present in supernovae (SNe) and PWNe which are supposed 
to be the sources of cosmic rays in the Galaxy. They are also 
present in relativistic flows ejected from ANGs. The astro- 
nomical phenomena are widespread and show synchrotron 
radiation produced by energetic particles. The accelerated 
particles have power-law spectra similar to cosmic rays. The 
energy of these shocks is the main energy of the system. 
There arc no other reasonable sources of cosmic rays but 
relativistic shocks. We claim that reflection process and 
Cracow acceleration produce almost all cosmic rays with 
energies higher than ~ 50 GeV. 



Appendix A: 

This Appendix will not be included in a review article. 

In 1997 I was carrying out simulations of relativistic 
shocks. I inspected the results of the simulations and no- 
ticed that the part of A9u caused by scatterings (A6 S ) 
diminishes with the time At u slower than the part of A6jj 
arising due to trajectory curvature in the uniform field com- 
ponent (A9i) when 7 grows and A is constant. I checked how 
quickly they are diminishing and got that A9 S decreases 
~ \/At u and A9i ~ At u . I understood at once that it re- 
sults from that the scattering is the diffusion process and 
the wandering in the uniform field is the linear process. I 
also known that the scattering enables the particles the re- 
turn to the shock and the wandering in the uniform field 
does not. 

I wanted to write it down in Bednarz & Ostrowski 1998, 
but M. Ostrowski did not want. He made up the general 
conclusion 'The phenomenon of decreasing a to at con- 
stant A and for growing 7 results from slower diminishing 
of the part of A9jj caused by scattering in comparison to 
A9jj arising due to trajectory curvature in the uniform field 
component.' 

I agreed because M. Ostrowski was resolute and there 
was no free space in the article. I thought that physicists 
know at once that the scattering is the diffusion process and 
the wandering in the uniform field is the linear process. 

During the 16th European Cosmic Ray Conference in 
Alcala de Hcnares in 1998 Y. Gallant and other man came 
up to me. Gallant wanted to know how the particles are 
accelerated in relativistic shocks. They did not understand 
the general conclusion in Bednarz & Ostrowski 1998. My 
English was very poor but I wrote down some mathematical 
formulae and an example and explained about the diffusion 
process, the linear process and what enables particles the 
return to the shock. Gallant got the point. He wrote the 
paper Gallant & Achterberg 1999. 
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